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ABSTRACT 

We discuss a model for the early assembly of supermassive black holes 
(SMBHs) at the center of galaxies that trace their hierarchical build-up far up in 
the dark halo 'merger tree'. Motivated by the observations of luminous quasars 
around redshift z ^ 6 with SMBH masses ~ 10^ Mq, we assess the possibility of 
an early phase of stable super-critical quasi-spherical accretion in the BHs hosted 
by metal free halos with virial temperature Tvir > lO^K. We assume that the first 
'seed' black holes formed with intermediate masses following the collapse of the 
first generation of stars, in mini-halos collapsing at z ~ 20 from high-cr density 
fluctuations. In high redshift halos with Tvir > lO^K, conditions exist for the 
formation of a fat disc of gas at Tgas ~ 5000 — lOOOOK. Cooling via hydrogen 
atomic lines is in fact effective in these comparatively massive halos. The cooling 
and collapse of an initially spherical configuration of gas leads to a rotationally 
supported disc at the center of the halo if baryons preserve their specific angular 
momentum during collapse. The conditions for the formation of the gas disc and 
accretion onto a central black holes out of this supply of gas are investigated, 
as well as the feedback of the emission onto the host and onto the intergalactic 
medium. We find that even a short phase of super-critical accretion eases the 
requirements set by the 2; ~ 6 quasars. 

Subject headings: cosmology: theory - black holes - galaxies: evolution - quasars: 
general 



1. Introduction 

The strongest constraint on the high-redshift evolution of supermassive black holes 
(SMBHs) comes from the observation of luminous quasars at z ^ 6 in the Sloan Digital Sky 
Survey (SDSS, Fan et al. 2001). The luminosities of these quasars, well in excess of 10^^ 
erg/s imply that SMBHs with masses ~ IO^Mq are already in place when the Universe is 
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only 1 Gyr old. The highest redshift quasar currently known, SDSS 1148+3251, ai z — 6.4, 

has estimates of the SMBH mass in the range (2 - 6) x IO^Mq (Barth et al. 2003, Wihott, 
McLure & Jarvis 2003). Among the seed BHs proposed, the less exotic (e.g. PopIII star 
remnants, gravitationally collapsed star clusters) are in the range 10^ — IO^Mq, forming at 
2; = 30 or less. 

To grow such seeds up to 10^ Mq requires an almost continuous accretion of gas, assuming 
Eddington accretion and a standard thin disc radiative efficiency for a Schwarzschild black 
hole, e ~ 0.1. If accretion is via a gaseous thin disc, though, the alignment of a SMBH with 
the angular momentum of the accretion disc tends to efficiently spin holes up (Volontcri 
et al 2004), and radiative efficiencies can therefore approach 30% (assuming a "standard" 
spin-efficiency conversion) . The accretion of mass at the Eddington rate causes the BH mass 
to increase in time as 



where t-^dd = 0.45 Gyr. Given M(0), the higher the efficiency, the longer it takes for the BH 
to grow in mass by (say) 20 e- foldings (Shapiro 2005). 

Yu & Tremaine (2002), Elvis, Risaliti, & Zamorani (2002), and Marconi et al. (2004), 
compared the local MBH mass density to the mass density accreted by luminous quasars, 
showing that quasars have a mass-to-energy conversion efficiency e ^0.1 (an argument 
originally proposed by Soltan 1982). This high average radiative efficiency, though, describes 
the population of SMBHs at redshift z < 5, and nothing is known about the radiative 
efficiency of pregalactic quasars in the early Universe. The numbers quoted above, in fact, 
may suggest that the picture that we have of the low redshift Universe may not apply at 
earlier times. 

In this paper we follow earlier papers in considering a scenario for the hierarchical 
assembly of SMBHs that traces their seeds back to the very first generation of stars, in 
mini-halos above the cosmological Jeans mass collapsing at 2; ~ 20 from the high-cr peaks of 
the primordial density field. 

However we introduce two new features: (i) we consider more explicitly the configuration 
of the gas from which the accretion occurs: a fat, dense disc of cold gas, likely to form in 
halos with Tyir > lO^K and zero metaUicity (Oh & Haiman 2002). And (ii) we explore the 
assumption that at early stages, accretion occurs at a 'Bondi' rate that is higher than the 
standard rate for e ~ 0.1 

The evolution of the first MBHs into the first pregalactic quasars and their impact on 
the re-ionization of the universe in this scenario has been explored recently by Madau et al. 
2004. We here expand and deepen the previous study, focusing in particular on the interplay 
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of gas cooling and BH feeding. 

In the next section we review the model for assembly of MBHs in cold dark matter 
(CDM) cosmogonies. We then address the conditions of the gas around the MBHs in the 
cores of high redshift mini-halos (§3), before and after the onset of pregalactic quasar activity. 
We compute the global evolution of the MBH population and discuss its imphcations in §4 
and finally summarize our results in §5. Unless otherwise stated, all results shown below 
refer to the currently favoured ACDM world model with Qm — 0.3, fl\ — 0.7, h — 0.7, 
fib = 0.045, as = 0.93, and n = 1. 

2. Assembly of pregalactic MBHs 

The main features of a plausible scenario for the hierarchical assembly, growth, and 

dynamics of MBHs in a ACDM cosmology have been discussed by Volonteri, Haardt, & 
Madau (2003), Volonteri et al. (2004), Madau et al. (2004). "Seed" holes with intermediate 
masses form as end-product of the very first generation of stars. They form in isolation within 
mini-halos above the cosmological Jeans mass collapsing at z ^ 24 from rare u-a peaks of 
the primordial density field (Madau & Rees 2001). Pregalactic seed IMBHs form within the 
mass ranges 20 < m. < 70 Mq and 130 < m. < 600 M© (Fryer, Woosley, & Heger 2002, 
Omukai & Palla 2003); for simphcity, within these intervals the differential black hole mass 
function is assumed to be flat. As our fiducial model we take u — A, corresponding to a mass 
density parameter in MBHs of Q, ~ 10~^. This is much less than the mass density of the 
supermassive variety found in the nuclei of most nearby galaxies, I^smbh = (2.1 ±0.3) x 10~^ 
(Yu & Tremaine 2002). Note that this choice of seed MBHs occupation is similar to that of 
Volonteri et al. 2003 (seed holes in 3.5a peaks aX. z = 20). 

3. Super-critical accretion and rapid growth of SMBHs 

3.1. Fat cold discs in high redshift haloes 

Metal-free gas in halos with virial temperatures Tvir > lO^K can cool in the absence of 
H2 via hydrogen atomic lines to ~ 8000 K. Oh & Haiman (2002) considered the cooling and 
collapse of an initially spherical configuration of gas in a typical halo with Tyir > lO^K. In the 
absence of H2 molecules, the gas can contract nearly isothermally at this temperature. Since 
coohng via the Lya hue has a very sharp cutoff at T < lO^K (Spitzer 1978) a 'thermostatic 
mechanism' operates: if the gas cools below ~ lO^K, the gas recombines and the cooling time 
rapidly increases. The gas then heats and contract until atoms are coUisionally re-ionized. 
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Halos, though, possess angular momentum, J, that can be related to the spin parameter 
A = J\E\^/'^/GMI'^, where E and Mh are the total energy and mass of the halo. The baryons 
probably preserve their specific angular momentum during collapse (Mo, Mao & White 1998), 
and settle into a rotationally supported disc at the center of the halo (Mo, Mao & White 
1998, Oh & Haiman 2002). 

Following Oh & Haiman (2002), we assume that a fraction fd of the gas settles into an 
isothermal, exponential disc, with gas temperature Tg^s, embedded in a dark matter halo 
of virial temperature Tvir, described by a Navarro, Frcnk & White (1997, hereafter NFW) 
density profile. The mass of the disc can therefore be expressed as Mdisc = fd{^b/^M)Mh. 

Under these assumptions, the disc scale radius is given by Rd ~ 2~^/^Ai?vir- 

For an isothermal exponential disc of radial scale length Rd, the number density of 
hydrogen at radius r and at vertical height z can be written as 



where no is the central density and zo is the vertical scale height of the disc at radius r. 
The central number density of the gas can be written as: 



Here Mh,g is the halo mass in units of 10^ M©, Tgas.sooo is the gas temperature in units of 
8000K, i?vir,6 is the virial radius in units of 6kpc, fd,o.5 is the fraction of gas in the disc 
normalized to 0.5, and Aq.os the spin parameter normalized to 0.05. Given the scaling of 
halo mass and virial radius in the spherical collapse model, no oc fd(Tgas/T^i^)~^ , so 
Zo oc A^/^^(Tgas/Tvir)^^^ and the disc is 'fat' when Tgas ~ Tvir, i-e. for haloes with virial 
temperature Tyir ^ lO^K whose gas can cool to ~ 8000 K due to hydrogen atomic line 
coohng. The disc thins out if the gas cools to lower temperatures in presence of metals or 
H2 molecules. 

Oh & Haiman showed that if only atomic line cooling operates, the discs arc stable to 
fragmentation in the majority of cases, and this fat disc contains essentially all the baryons. 
In the rare cases where fragmentation is possible, there would presumably then be star 
formation. However it seems inevitable that a substantial fraction of the baryons will remain 
in the fat disc at 5000K < Tgas < lO^K (at the bottom of this range if there were no significant 
heat input, but shifting towards the peak of the HI cooling curve if there is substantial 
heating) . 






-3 



(3) 



- 5 - 



3.2. Accretion onto a central MBH 

As described in §2, we assume that one seed MBH forms in each of the rare density 
peaks above 4-cr at z = 24. The MBH is supposed to be at rest in the center of the halo, 
after the formation and collapse of a PopIII star in the central region of the mini-halos. 

Let us now consider a MBH seed in a metal free halo with Tvir ^ 10"^ K. As delineated 
in the previous section, a fat disc of cold (5000K < Tgas < lO'^K) gas likely condenses at the 
center of the halo. This gas can supply fuel for accretion onto a MBH within it. We will 
now discuss how stable super-critical accretion (i.e. much larger than the Eddington rate) 
of the gas can proceed. 

In a spherical geometry, matter can infall onto the hole at a rate that can be estimated 
using the Bondi-Hoyle formula (Bondi & Hoyle 1944): 

-'WBondi — 5 14J 

where a is a dimensionless parameter of order unity and Cg ~ 10(T/10^K)^/^ kms"^ is the 
sound speed of the gas. 

It is useful to compare the above accretion rate to the Eddington rate: 

^Bondi _ ^^^^ y-3/2 



MEdd 

where the central density no,4 is in units of lO^cm"^ and the MBH mass, Mbh,3, is normalized 
to 10^ Mq. Clearly for the typical black hole masses and gas densities we are considering 
here, Meondi >> ^Edd- 

The Bondi radius for the MBHisra^^c = GMbh/c^ ^ 0.05 pc(Mbh/10^ Mo)(rgas/8000K)-^ 
Let us now compare the Bondi accretion radius to the fat disc thickness: 



6 X 10"^ Mbh,3 ^gasjooo ^0,1 Km- 



For MBHs in halos with Tyir ^ 10^ K, the accretion radius is then well within the disc, 
surrounded by gas with almost constant density uq and thus a quasi-spherical geometry. In 
particular, if the gas disc rotates as a rigid body, the transverse velocity at the Bondi radius 
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is very small compared to the free fall velocity. The latter is, by definition, equal to the 
sound speed, while the former is a factor rg^/Rci smaller than vd-, the rotation velocity of 
the gas disc at i?^, vd ^ Cg. The angular momentum is nonetheless too large for the gas to 
fall radially into the hole and a tiny accretion disc forms. The outer edge of the accretion 
disc, Tin, can be expressed as: 

^ = 3 X 10^ vl^,, Ml^^, R;1, (6) 

where rs is the Schwarzschild radius of the hole, vr, is in units of 10 kms~^ and we have as- 
sumed specific angular momentum conservation inside race- The accretion disc has therefore 
a size of order of the trapping radius: 



rtr = rs^l — = 0.12 r,^ Vj^^q M^^ 3 i?^ir,6^o,4, (7) 

i.e. the radius at which radiation is trapped as the infall speed of the gas is larger than 
the diffusion speed of the radiation. The trapping radius is equivalent to the radius of 
spherization defined by Shakura & Sunyaev (1973), where the thickness of the disc becomes 
of the same order as rtr; at smaller radii the the inflow is quasi-spherical. 

Begelman (1979) and Begelman & Meier (1982) studied super-critical accretion onto a 
BH in spherical geometry and quiescent thick discs respectively. In the spherical case, radia- 
tion pressure cannot prevent the accretion rate from being super-critical, while the emergent 
luminosity is limited to ^ -^Edd- The case with angular momentum is more complicated 
(especially with regard to the role of winds). Though this issue remains unclear, it still 
seems possible that when the inflow rate is super-critical, the radiative efficiency drops so 
that the hole can accept the material without greatly exceeding the Eddington luminosity. 
The efficiency could be low either because most radiation is trapped and advected inward, or 
because the flow adjusts so that the material can plunge in from an orbit with small binding 
energy (Abramowicz & Lasota 1980). 

Despite the uncertainties it seems worthwhile to explore the consequences for black hole 
growth of an early phase when the MBH can accept most the mass infalling at a rate that 
can be estimated using the Bondi-Hoyle formula, and can therefore grow much more rapidly 
than the Eddington rate would allow. 

The conditions for quasi-spherical geometry hold when the accretion radius is much 
smaller than the disc thickness, i.e. when the MBH and the halo are small and the gas 
temperature is ~ 8000 K. The gas temperature will not be much larger than ~ 10^ K, as the 
dense gas can radiate away all the energy injected by the quasar in photoionizing photons. 
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On the other hand, if the gas is able to cool down to lOK < Tgas < 200K, due to 
formation of molecular hydrogen or metal pollution, collapse and fragmentation in the disc 
will suddenly halt the Bondi-style accretion. We therefore argue that this super-critical 
accretion phase would end when the Universe is enriched by metals at 6 < ^ < 10. 

In the absence of metals, H2 would be the main coolant, but the strong internal feedback 
provided by the pregalactic quasars emission within a galaxy easily suppresses H2 formation. 

In fact, assuming a nonthermal power-law component (a = 1), in the 11.15-13.6 eV range, 
the total rate of dissociations produced by the pregalactic quasars in the gas disc is A'diss ~ 
10^^ (^bh/10'^ M0)(no/lO'^cm~^)(2;o/pc) s~^. The total number of molecules in the disc is 
A^Ha ^ 10®^ (Miisc/10^ M©). The timcscale for dissociating all the H2 molecules present in 
the disc is thus very short. While the pregalactic quasar is shining, H2 formation is therefore 
completely halted in the disc, and cooling below ~ 4000K is impossible. 



4. Results 

To assess the relative importance of this processes, we trace the evolution of a MBH 
along with its host halo. We generate Monte Carlo realizations (based on the extended 
Press-Schechter formalism) of the merger hierarchy of a — 10^^ Mq halo at 2; = 6. The 
halo mass is chosen by requiring that the number density in haloes more massive than 
matches the space density of quasars at 2; = 6 (Fan et al. 2003, 2004). We then extract from 
the trees the mass-growth history of the main halo and of smaller satellites (selecting only 
density peaks above 4-(7) and grant these halos a MBH seed. 

When the halo virial temperature becomes larger than lO^K, the MBH is surrounded 
by a coldish (Tgas = 8000K) disc. The accretion rate is set at the Bondi rate, as determined 
by the gas density at the accretion radius. The luminosity, though, does not exceed the 
Eddington luminosity. During the accretion process, the radius of the inner disc increases 
steeply with the hole mass, thus making super- Eddington accretion less likely to be sustained. 
We assume here that, when the radius of the accretion disc becomes a factor of 5 larger 
(though the choice of the exact value is somehow arbitrary) an outflow develops, blowing 
away the disc. A subsequent major merger then has to be awaited before a fresh supply 
of fuel is available to the MBH. Typically, the conditions for super-critical accretion are 
satisfied only once per MBH. If they are not, the MBH accretes at the Eddington value 
during the susequent accretion episode(s). The radiative efficiency then evolves with the 
MBH spin, adopting the standard definition for circular equatorial orbits around a Kerr 
hole. The MBHs spin is modified during the accretion phase as described in Volonteri et al. 
2004. Figure 1 shows the evolution of MBH (mass and accretion rate) and disc in the main 
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halo of one of the merger trees. This halo represents a S-cr density fluctuation at 2; = 24. 
The virial temperature of the main halo is larger than lO^K at ^ ~ 24 and the MBH starts 
growing very early The accretion rate is initially super-critical by a factor of 10 and grows 
up to a factor of about 10*^, thus making the flow more and more spherical (see Equation 
7). On the other hand, the whole 'plump' accretion disc grows in size until it crosses the 
trapping surface and reaches the assumed threshold for the end of the activity. The figure 
shows also the effect of a larger or smaller fraction of the baryons ending up in the fat disc. 
A large supply of fuel (/^ = 0.5) triggers rapid accretion early on, but also the accretion is 
quenched sooner. In both cases the spin parameter is set as A = 0.05. MBHs in halos with 
small spin parameter have an early and short super-critical accretion episode, (cfr. eq. 3 
and 6) while in case of high spin parameter, the super-critical phase happens, typically with 
a more dramatic growth of the BH mass, at lower redshift. 

Our assumed threshold for MBH formation (density peaks above A-a) selects only the 
most massive and rare halos at a given redshift. This ensures that super-accreting systems 
are not widespread. Let us consider comparatively smaller halos (e.g. massive satellites of 
the main halo. Figure 2), with virial temperature below IC^K at z ~ 24. They accrete mass 
along the merger hierarchy, until their mass implies a virial temperature larger than 10"^K. 
Atomic cooling, therefore, becomes effective at lower redshift, so that there is only a short 
time (if any) for the rapid growth to take place. The effect of supercriticality on black holes 
hosted in smaller halos is thus mild. Note that we have created a whole set of merger trees, 
and they give qualitatively very similar results. 



5. Discussion 

In this paper, we have envisaged an early stage of super-critical accretion during the 
global evolution of a SMBH (see also Kauffmann & Haehnelt 2000). Fuel is supplied by 
a dense gaseous disc forming in halos with Tyi^ > lO^K, where hydrogen atomic coohng is 
effective. 

If the disc rotates as a rigid body, the gas infalling on the central MBH has small 
transverse velocity. A tiny accretion disc can form within the radius at which radiation is 
trapped and the MBH can accept most of the infalling mass. Even if this phase lasts only 
until the universe in enriched by metals, the quick start allows the holes in the most massive 
halos to reach the high SMBH masses suggested by the SDSS quasars. It is worth noting that 
super-critical accretion in the fashion described here occurs at very high redshift, and that 
stops well before 2; ~ 6. We do not expect, therefore, that SDSS quasars are accreting above 
the Eddington value, consistently with observations (Barth et al. 2003, Willott, McLure & 
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Jarvis 2003). A signature of super-critical accretion is the occurrence of outflows, that we 
envisage would quench activity. Outflows from these sources could leave their imprint by 
spreading metals into the IGM early-on. 

One might be concerned that this model overproduces large SMBHs at low redshift 
This is not the case, since only a tiny fraction of halos have Tyir ^ 10^ K, when the Uni- 
verse is still metal free. Assuming that metal pollution starts affecting the disc cooling 
and fragmentation at z < 10, only a fraction ^ 2 x 10~^ of halos at this time is likely to 
be the outcome of a merger hierarchy which involves a seed MBH aX z > 20. So, even if 
all these halos with mass Mh ~ lO^'' M© hosted a MBH with mass ~ 10^ M0, the density 
i^MBH = (0.002 flM {Mbb)) / Mh = 10~* would still be much lower than the local one. We will 
discuss this issue in detail, along with the possibihty that black holes can be displaced from 
galaxy centers and ejected into the IGM by the 'gravitational rocket' effect in a subsequent 
paper. 

As a side result, this model allows a very 'economic' re-ionization in terms of BH seeds 
density. Madau et al. (2004) found that pregalactic quasars powered by MBHs forming 
in 3.5(7 peaks will reionize the IGM if they accrete at the Eddington rate a mass of order 
10~^ the mass of the host halo in every major merger. Alternatively, seed holes must be 

more numerous (e.g. 3-0" density fluctuations) at the start in order to sustain the early 
production of ionizing radiation. This latter assumption would imply that all galaxies in 
the local Universe with total mass > 5 x IO'^Mq are expected to host a central MBH. 
Allowing for a super-critical accretion phase, the number of ionizing photons per hydrogen 
atom produced by pregalactic quasars approaches unity at 2; > 15, even assuming a very low 
initial density in seeds (e.g. 4-0" peaks), thus relieving the need for a very large number of 
seeds - difficult to reconcile with low redshift constraints. 

The growth of SMBHs that can power SDSS quasars can be explained within a ACDM 
universe also assuming a more optimistic view in terms of accretion and merging. Yoo 
& Miralda-Escude (2004) showed 2; ~ 6 quasars can be explained assuming continued 
Eddington-hmited accretion onto MBHs forming in halos with Ty^^ > 2000K at 2; < 40. 
Their model assumes, also, a much higher influence of BH mergers in increasing the MBH 
mass: a contribution by itself of order 10^ M©. Their investigation takes into account the 
negative feedback that dynamical processes at BH mergers ('gravitational rocket', see also 
Haiman 2004) impose on SMBHs growth. A heavy influence of MBH mergers in building-up 
SMBHs can be probed by their gravitational waves emission, by the planned Laser Interfer- 
ometer Space Antenna (LISA). 

As pointed out by Shapiro (2005), explaining the presence of SMBHs at 2; ~ 6 is 
extremely difficult if accretion occurs at the Eddington rate and via a standard thin disc. If 
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alignment between the accretion disc and the MBH spin is efficient on very short timescales, 
MBHs quickly end-up maximally spinning, with a radiative efficiency e ~ 0.4 (cfr. Equation 
1). Shapiro suggests therefore that accretion occurs via discs where viscosity is due to 
magneto-hydrodynamical (MHD) turbulence. Gammie et al. (2004) simulations suggest in 
fact that the maximum spin MBHs can achieve by coupling with MHD discs is smaller than 
1 and the corresponding maximum radiative efficiency is e ~ 0.19. The gain in time available 
for SMBH growth is less than a factor of 3, though, thus probably requiring that dynamical 
negative feedback is not important. 

The issues we have discussed in this paper will be clarified by better evidence on mini- 
quasars and on the re-ionization history at 2; ~ 10. Future missions such as LOw Frequency 
ARray (LOFAR) and the Square Kilometer Array (SKA) will probe directly the IGM up to 
high redshift. Our discussion also provides added motivation for ongoing studies of the flow 
patterns that occur when gas falls towards a hole at a supercritical rate. 
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Fig. 1. — Evolution of the accretion rate and MBH mass growth in the main halo of the 
merger trees, = 10^^ Mq at z = 6. Top: fd = 0.5. Bottom: = 0.1. Left panel: 
MBH mass as a function of redshift for the model discussed in this paper {solid line) and 
assuming Eddington accretion rate at all times, with e = 0.15 {dashed line). Right panel: 
MBH accretion rate, in units of the Eddington rate. 
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Fig. 2. — Same as Figure 1, but for a massive satellite halo. This halo is a 4-0" peak in 
the density fluctuations fleld, and represents therefore the minimum mass of a halo which is 
granted a seed MBH in our model. 



